EPJ manuscript No. 

(will be inserted by the editor) 



Jet Tomography in the Forward Direction at RHIC 

G.G. Barnafoldii2 p, Levai\ G. Papp^, and G. Fai^ 



O 
O 

(N 

Oh' 
(D 

in 

(N 

(N 
> 
m 
(N 
O 

a^ 
o 
\o 
o 

Oh! 



^ RMKI KFKI, P.O. Box 49, Budapest 1595, Hungary 
^ Dept. of Physics of Complex Systems, Eotvos University, Pazmany P, 



X 



1/A, Budapest 1117, Hungary 
Dept. for Theoretical Physics, Eotvos University, Pazmany P. 1/A, Budapest 1117, Hungary 
Center for Nuclear Research, Kent State University, Kent, OH 44242, USA 

Received: date / Revised version: date 

Abstract. Hadron production at high-pT displays a strong suppression pattern in a wide rapidity region 
in heavy ion collisions at RHIC energies. This finding indicates the presence of strong final state effects for 
both transversally and longitudinally traveling partons, namely induced energy loss. We have developed a 
perturbative QCD based model to describe hadron production in pp collision, which can be combined with 
the Glauber - Gribov model to describe hadron production in heavy ion collisions. Investigating AuAu and 
CuCu collisions at energy ^/s — 200 AGeV at mid-rapidity, we find the opacity of the strongly interacting 
hot matter to be proportional to the participant nucleon number. Considering forward rapidities, the 
suppression pattern indicates the formation of a longitudinally contracted dense deconfined zone in central 
heavy ion collisions. We determine parameters for the initial geometry from the existing data. 

PACS. 12.38.Mh Quark-gluon plasma in quantum chromodynamics - 24.85+p Quarks, gluons, and QCD 
in nuclei and nuclear processes - 25.75-q Relativistic heavy ion collisions 



1 Introduction 

One of the most interesting challenges in high-energy heavy 
ion collisions is to extract the properties of the produced 
hot dense medium. The discovery of strong pion suppres- 
sion at RHIC energies in AuAu collisions ^12113] and its 
detailed analysis was the first milestone along the 

way to determine the parameters of the central space- 
time region of heavy ion collision. Applying the so called 
"jet-tomography" method [HltZIISIini one can "measure" the 
averaged density of the color charges in the deconfined re- 
gion on the basis of the induced energy loss of partons 
traversing the dense region. Previously we have investi- 
gated the induced energy loss for transversally propagat- 
ing partons 10, analyzing mid-rapidity data at different 
centralities in AuAu collisions JT,JJj and we have deter- 
mined the opacity parameter of the hot dense matter in 
the transverse direction. The next level is to perform to- 
mography into the forward (or backward) directions and 
to extract integrated information about the 3-dimensional 
structure of the color charge density. 

In this paper we display forward tomography in AuAu 
collisions at RHIC energy. Our aim is to understand the 
latest BRAHMS result ^31 obtained in AuAu collisions, 
where the measured nuclear modification factor, RAAipr), 
displays similar values at different rapidities (y = and 
y = 3.2). The explanation of this requires the understand- 
ing of competing nuclear effects, namely shadowing, mul- 
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tiscattering, induced energy loss, and the qualitative con- 
sideration of the geometry of the dense partonic region. 
On the other hand, these data offer a unique possibility 
to extract information about the properties of the dense 
matter in the longitudinal direction in the very early stage 
of the heavy ion collisions. 

Our survey will start with the pp and the dAu data 
as baseline. At large forward rapidities the indication of 
increasing suppression has been seen in dAu collisions in 
pion and charged hadron production [T^ItT^ . The analysis 
of these data helped us include the appropriate shadowing 
at high rapidities 16 ,17 . Thus we are able to determine 
quantitatively the appropriate energy loss in AuAu colli- 
sions, assuming that all final state effects are connected 
to induced radiative energy loss and the geometry of the 
dense matter. 

We briefly introduce our perturbative QCD improved 
parton model [TnilTHj and the description of different nu- 
clear effects. The initial state nuclear shadowing and mul- 
tiple scattering effects are able to reproduce the tt" spec- 
trum in dAu collisions |16II19| . The final state radiative 
parton energy loss is modeled by the GLV- method . We 
extract the opacity of the dense matter in different cen- 
trality regions in AuAu and CuCu collisions at ^/s ~ 200 
AGeV at mid-rapidity. On the basis of the extracted opac- 
ity values we perform the forward tomography calculations 
and investigate existing AuAu data at large forward ra- 
pidities. 
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We discuss the obtained results and connect them to 
the longitudinal geometry of the dense region. Our anal- 
ysis is based on the application of radiative energy loss of 
jets in the mid-rapidity and forward rapidities. The role of 
other nuclear effects (e.g. possible appearance of scatter- 
ing energy loss for light quarks PUII^ I will be discussed 
in a forthcoming paper. 



2 Theoretical Description of the Model 

The high-pT inclusive pion spectrum in a heavy ion colli- 
sion can be calculated in a pQCD-improvcd parton model. 
Originally this was developed for nuclcon-nucleon colli- 
sion (basically pp), and extended by a Glauber- type colli- 
sion geometry and initial state nuclear effects for nucleus- 
nucleus, A A' coUisions as P^23n24| : 
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Here tA{b) = J dzpA{b,z) denotes the nuclear thickness 
function, and it is normalized as usual: J d btA{b) = A. 
For proton and deuteron we use a sharp sphere approx- 
imation, while for heavy nuclei the Wood-Saxon formula 
is applied for the nuclear density distribution, pA{b,z). 

In our next-to-leading order (NLO) calculation, da/dv 
represents the Born cross section of the partonic sub- 
process and Kab,cis,v,w,Q,QR,QF) is the correspond- 
ing higher order correction term, see Ref.s |22II23II24| . We 
fix the factorization and renormalization scales and con- 
nect them to the momentum of the intermediate jet, Q = 
Qr = (4/3)pq (where pg = pt/zc), reproducing pp data 
with high precision at high pT jl9| . 

To take into account the transverse momentum dis- 
tribution, we defined the following 3-dimensional parton 
distribution function (PDF): 

fa/p{Xa,kTa,Q^) = fa/p{Xa,Q^) ■ Oa/pO^Ta) ■ (2) 

Here, the function fa/p{xa,Q^) represents the standard 
longitudinal NLO PDF as a function of momentum frac- 
tion of the incoming parton, Xa, at the factorization scale 
Q. In the present calculation we use the MRST(cg) pa- 
rameterization |25| . The partonic transverse- momentum 
distribution in 2 dimensions, gg^/plkx), is characterized by 
an "intrinsic fc^" parameter as in Refs. |18ll24j . In our phe- 
nomenological framework we assumed a Gaussian func- 
tion [IHIEEI- 



Nuclear multiscattering is accounted for through a broad- 
ening of the incoming parton's transverse momentum dis- 
tribution function, namely an increase in the width of the 
Gaussian: 



{k^)pA — {k'^)pp + C ■ hpA{b) 



(3) 



Here, {k'^)pp = 2.5 GeV'^ is the width of the transverse 
momentum distribution of partons in pp collisions |18lll9j , 
hpA{b) describes the number of effective NN collisions at 
impact parameter 6, which impart an average transverse 
momentum squared C. The effectivity function hpA{b) can 
be written in terms of the number of collisions suffered by 
the incoming proton in the target nucleus. In Ref. |18| 
we have found a limited number of semihard collisions, 
3 < J^m < 4 and the value C = 0.35 GeV^. 

We take into account the isospin asymmetry by using 
a linear combination of p and n PDFs. The applied PDFs 
are also modified inside nuclei by the "shadowing" effect, 
applying the parametrization of Ref. '27' . 

The last term in the convolution of eq. is the frag- 
mentation function (FF), D^{zc, Q^)- This gives the prob- 
ability for parton c to fragment into a pion with momen- 
tum fraction Zc at fragmentation scale Q — (4/3)pT- We 
apply the KKP parameterization |28j . 

Finally we include jet-quenching as a final state effect. 
The energy loss of high-energy quarks and gluons traveling 
through dense colored matter can measure the integrated 
density of the colored particles. This non-Abelian radia- 
tive energy loss, AE{E, L), can be described as a function 
of parton opacity, (or mean number of jet scatterings), 
n = L/A, where L is the interaction length of the jet and 
A is the mean free path in non-Abelian dense matter. In 
" thin plasma" approximation the energy loss to first order 
is given by the following expression |S]: 
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Here Cr denotes the color Casimir of the quark or gluon 
jet, p/ X ~ aj.ppart is the transport coefficient of the medium, 
which is proportional to the parton density, Ppart, and p 
denotes the color Debye screening scale. N{E) is an energy 
dependent factor with asymptotic value 4. 

Considering a time-averaged, static plasma, the aver- 
age energy loss, AE, will modify the argument of the FFs: 



zl D^/c{z*,t 
Zr 7rz2 



(5) 



Here z* = Zc/ (1 — AE/pc) is the modified momentum 
fraction. 
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Fig. 1. The nuclear modification factor for vr'' production in 
AuAu collision, RauAu{pt), at energy ^/s = 200 ^GeV at mid- 
rapidity, measured by PHENIX 29 30 . The solid lines indicate 
our best fit with the applied opacity parameters, L/X. 



We will present results on hadron productions through 
the nuclear modification factor, RAA'{pT,b) defined by 



RAA'{PT,b) 



1 



(6) 



Here is the average number of binary collisions, as a 
function of px at different impact parameter ranges. 



3 Jet Tomography in Different AA Systems 

The opacity parameter, n = L/X, can be determined by 
finding the best fit for energy loss and comparing the the- 
oretical results to the data points on the nuclear modifica- 
tion factor. In this section we present our jet tomography 
results in the mid-rapidity region. 

Fig. n displays our fit on opacity in AuAu collision 
at different centralities, using the preliminary data of the 
PHENIX collaboration on tt" production at mid-rapidity 
at ^ = 200 AGeV energy PMSU] . Experimental data in- 
dicate a stronger suppression pattern at > 4 GeV/c, 
thus we read the L/X values for this region. Here the nu- 
clear modification factor is flat, as Fig. ^ displays. The 
extracted opacity values are boxed at the top of the pan- 
els, together with the centralities. The plotted RAuAuiPT) 
are calculated at the mean values of the opacity, the errors 
of the fits are approximately ±0.25. 

The measured suppression pattern in AuAu collision is 
reproduced well in all centrality bins. In the most central 
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Fig. 2. The nuclear modification factor for pion produc- 
tion in CuCu collision, RcuCu{Pt) , in mid-rapidity at energy 
^ = 200 ^GeV. Data are from PHENIX (7r°) ED and STAR 
( '' +'^~ ) 1^ . The solid lines indicate our best fit with the 
applied opacity parameters, L/X. 



(0 — 10%) collisions we have obtained a maximal opacity 
value, L/X — 4.0±0.25. In the most peripheral collisions 
(80-92%) we have extracted the opacity L/X = 1.0±0.25. 

We analyzed the data from CuCu collisions at ^/s = 
200 ^GeV energy in a similar fashion. Fig. |21 displays the 
existing data and our calculations. We have found prelimi- 
nary data on tt^ production from the PHENIX Collabora- 
tion in the 0-10% and 10-20% centrality bins [HI- We are 
able to calculate pion spectra with high precision in our 
pQCD framework, thus we compare data and theory di- 
rectly in these two centrality bins. In parallel, STAR data 
exist on unidentified charged hadrons [K^ := ^ 1^^ ) for 
more centrality bins |32| . In the most central bins we see 
that the nuclear modification factors overlap for identified 
pions and unidentified charge hadrons in the very high-p^ 
region, pr > 6 — 8 GeV/c. (This indicates the restora- 
tion of the fragmentation limit for proton and antipro- 
ton production in this region, thus the validity of pQCD 
based calculation for all hadron species — for smaller mo- 
menta the contribution from quark coalescence should be 
included for baryon and antibaryon production 33 3^). 
Thus, focusing on the highest momentum region, we ex- 
tracted the opacity values from the STAR data on uniden- 
tified charge hadrons and display the data and the theoret- 
ical results on Fig. El for two more centrality bins, namely 
for 20 - 30% and 30 - 40%. 

Fig. 13 summarizes all mid-rapidity results on opacity 
for AuAu collisions in the — 92% centrality region {tri- 
angles) and for CuCu collisions in the — 40% region 
(squares). The opacity parameter, which is connected to 
the color charge density, is decreasing approximately lin- 
early with increasing centrality, and the slopes of the opac- 
ity curves are very similar for AuAu and CuCu collisions. 
Furthermore, the mid-central (30 — 40%) AuAu collisions 
are characterized by the same opacity parameter as the 
most central (0 — 10%) CuCu case. This overlap can be 
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Fig. 3. The centrality dependence of the opacity parameter, 
L/A, in AuAu (triangles) and CuCu [squares) colhsions at 
energy ^ = 200 ylGeV. 
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Fig. 4. The opacity parameter, L/A, as a function of the 
participant number, (Npart) in AuAu (triangles) and CuCu 
(squares) colhsions at ^/s = 200 AGeV. 



seen clearly in the measured nuclear modification factors 
for AuAu and CuCu collisions at different centralities and 
it has been discussed in Refs. |35II36II37| . This effect may 
present a real challenge to theoretical models as well as 
the consequences of any explanation. 

Fig. 21 displays the extracted opacity parameters as a 
function of the participant number, (Npart)-, in AuAu col- 
lisions (triangles) and CuCu collisions (squares). Here the 
overlap of the two centrality dependent opacity parameter 
sets is clearly seen. Fitting the opacity values with a con- 
tinuous function of the participant number (see solid line 
on Fig. 2)1, the following functional form has been found: 

L/A = (0.62 ± 0.09) • (A^paw)°-''^°-°' ■ (7) 

This fit agrees well with the simple expectation based on 
geometry, namely L cc A^^^ cc N^^^'^f [55] . 

Assuming constant mean free path (A) and zero for- 
mation time, the above fit validates the idea of a volume- 
dominated jet energy loss for the whole volume of the hot 
dense matter. This was assumed in the GLV-method 6 , 
thus our analysis is self-consistent. Combining this result 
with the expression of the induced energy loss in eq. , 

one can obtain the simple connection, e = AE/E oc N^l^^. 
In the forthcoming sessions we will assume the presence 
of volume-dominated energy losses. 

On the other hand, in models which lead to surface- 
dominated jet-energy loss j39ll4()ll41j or which assume non- 
zero formation times |42| the above picture is modified 
drastically. It would be interesting to see the centrality 
dependence in a repeated analysis with these assumptions 
and the appropriate plots equivalent to Fig. |3| and Fig. ^ 
However, this is not our aim, here, in this paper. 



4 Jet Tomography at Large Rapidities 

The jet-quenching pattern of the nuclear modification fac- 
tor is driven by the space and time averaged value of the 
color-charge density. We are interested in the presence of 
possible asymmetry between transverse and longitudinal 
directions seen by the particles produced in different ra- 
pidities. Since the produced hot dense matter is charac- 
terized by different dynamical behavior in transverse and 
longitudinal direction, we expect to see this difference dur- 
ing the analysis of available data. Here we analyze hadron 
production in the forward rapidity region using our theo- 
retical framework and the latest data at RHIC. 

In the latest BRAHMS results obtained in the most 
central AuAu collisions |13II14| the nuclear modification 
factor, Raa{pt), displays similar values at different ra- 
pidities, namely at y = and y = 3.2. We expect, that 
the combination of appropriate geometry of the hot dense 
partonic region and the competing nuclear effects, namely 
shadowing, multiscattering, and induced energy loss ex- 
plain this result. Thus, if we extend our pQCD based cal- 
culation and include these latter nuclear effects (which 
seems to be under control at higher rapidities [1^1117) 1. 
then the analysis of high rapidity data offers a unique 
possibility to extract information about the properties of 
the dense matter in the longitudinal direction in the very 
early stage of the heavy ion collisions. 

The BRAHMS data pHtTl] are available in the inter- 
mediate-pT region (pT < 4 GeV/c). We will assume that 
the suppression factors at different rapidities are approxi- 
mately the same in the high-p^ region and we will perform 
our calculations in this manner. This means a suppression 
factor ~ 5 for AuAu collisions at ^/s = 200 AGeV (seen in 
the most central collisions in the mid-rapidity [^) and a 
suppression factor 3 at = 62.4 AGeV (see PHENIX 
data in the most central AuAu collisions '30 .) 
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and CuCu collisions at energies ^/s — 200 AGeV. 



Fig. displays our results in the (pseudo)rapidity re- 
gions ry = 0.0, 1.0, 2.2 and 3.2 for central AuAu collisions 
at energies = 200 AGeV {left panels) and 6 2.4 AGeV 
{right panels). Data are from BRAHMS 13,14 at higher 
rapidities. The mid-rapidity data are from PHENIX |21 
1201 on pion and from STAR on charged hadron pro- 
duction. We display the extracted L/X opacity values for 
all rapidities. The extended bands indicate the theoretical 
uncertainties on the induced energy loss denoted by the 
errors in the opacity values. 

Fig. El summarizes the obtained opacity parameters in 
wide rapidity region for the most central (0 — 10%) AuAu 
collisions at ^/s — 200 AGeV and 62.4 AGeV energies. 
(We have included the opacity parameter from CuCu col- 
lisions at y/s = 200 AGeV.) In the highest energy domain 
the opacity parameter is decreasing with increasing ra- 
pidity from the maximal value of L/X = 4.0 ± 0.25 at 
mid-rapidity to L/X = 1.5 ± 0.25 at 77 = 3.2. Decreasing 
the collision energy to y/s = 62.4 AGeV a smaller opac- 
ity, L/X = 3.25 ± 0.25 has been extracted at 77 = 1. Since 
no data available at higher rapidities at this energy, we 
assumed the unmodified suppression factor. This way we 
have obtained L/X = 0.75 ±0.25 at 77 = 3.2. Fig.[Slcan in- 
dicate the validity of this assumption. Fig.Elclearly shows 
an approximately linear decreasing tendency for the opac- 
ity parameter as 77 is increasing. We have found that at 
large forward rapidities the interplay between a stronger 
shadowing and a linearly (in 77) decreasing quenching ef- 
fect is able to maintain a rapidity independent nuclear 
modification factor. 

Fig. ini displays one more interesting message, what we 
were looking for. Namely, comparing the extracted opacity 
values in the mid-rapidity and in the most forward rapid- 
ity, one can see a factor of 3 difference at both energies, 
Vs = 62.4 AGeV and ^ = 200 AGeV. According to our 
assumption the produced hot dense deconfined matter is 



homogeneous, thus we have a uniform A value (which may 
depend on energy but not on geometry). Thus the above 
result indicates that longitudinally traveling partons see 
less colored matter than those traveling in the transverse 
direction. This result can be explained on the basis of pure 
geometry combined with dynamics. 

Fig. displays a schematic picture for the time evolu- 
tion of the formed hot dense matter and the outcoming 
jets. Let us characterized the hot matter with the mean 
free path of A = 1 fm. A jet, created in the central re- 
gion of the collision and producing mostly mid-rapidity 
hadrons, is traveling transversally through an Lt ~ 4 fm 
length. This jet looses a large portion of its energy and 
indicates an opacity Lt/X — 4. Hadrons in the forward 
rapidities are produced from forward jets. These jets are 
moving mostly longitudinally and after passing a (con- 
tracted) thin region of the compressed matter (character- 
ized by an effective length oi L^ ~ 1.5 fm) they reach very 
quickly the longitudinally expanding surface. Both jet and 
expansion surface are moving with speed of light, thus the 
comoving jet does not loose more energy and an opacity 
Ll/X ~ 1.5 can be extracted at the highest rapidities. 



Am + Am 




t(fm) 



Fig. 7. Time evolution of the quark or gluon jet traveling 
through the expanding hot dense matter into the transverse 
and longitudinal direction. 
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Decreasing the collision energy, the density of the pro- 
duced hot matter is decreasing and the mean free path is 
increasing. We do not expect strong modification in the 
geometry and the comoving dynamics. Thus our results 
at i/i = 62.4 AGeV indicate a 25 % increase in the mean 
free path comparing to the case of i/s = 200 AGeV. 



5 Summary 

We have analyzed the measured suppression pattern in 
AuAu and CuCu collisions at energy ^/s — 200 AGeV at 
mid-rapidity. We have found the opacity of the strongly 
interacting hot matter to be proportional to the partici- 
pant nucleon number. This result indicates the presence 
of drastical microscopical mechanisms for energy stopping, 
which are working in a uniform way at given energy and 
producing similar energy densities in heavy ion collisions 
of different nuclei. 

Investigating hadron production in forward rapidities, 
the suppression pattern indicated the formation of a lon- 
gitudinally contracted dense deconfined zone in the most 
central heavy ion collisions. We could have determined the 
initial geometry from the existing data. The decrease of 
the opacity into the longitudinal direction is explained by 
the appearance of comoving dynamics. 

Our one-hadron tomography can be combined with di- 
hadron correlation studies 43 , 44 to obtain more detailed 
information about the geometry of the hot region. Energy 
and centrality dependence in CuCu collisions could be 
used to verify our results. 
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